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Abstract 

An axisymmctric model for the Crab Nebula is con- 
structed to examine the flow dynamics in the nebula. 
The model is based on that of Kennel and Coro- 
niti (1984), although we assume that the kinetic- 
energy-dominant wind is confined in an equatorial 
region. The evolution of the distribution function of 
the electron-positron plasma flowing out in the neb- 
ula is calculated. Given viewing angles, we reproduce 
an image of the nebula and compare it with Chandra 
observation. 

The reproduced image is not a ring-like but rather 
'lip-shaped'. It is found that the assumption of 
toroidal field does not reproduce the Chandra image. 
We must assume that there is disordered magnetic 
field with an amplitude as large as the mean toroidal 
field. In addition, the brightness contrast between 
the front and back sides of the ring cannot be repro- 
duced if we assume that the magnetization parame- 
ter a is as small as ~ 10~^. The brightness profile 
along the semi-major axis of the torus is also exam- 
ined. The non-dissipative, ideal-MHD approximation 
in the nebula appears to break down. 

We speculate that if the magnetic energy is released 
by some process that produce turbulent field in the 
nebula flow and causes heating and acceleration, e.g. 
by magnetic reconnection, then the present difficul- 
ties may be resolved (i.e. we can reproduce a ring 



image, and a higher brightness contrast). Thus, the 
magnetization parameter cr can be larger than previ- 
ously expected. 

stars: pulsars: general - ISM: individual: Crab 
Nebula 

1 Introduction 

A standard picture of the Crab Nebula was given by 
Kennel and Coroniti (KC; 1984). According to their 
picture, a super-fast magnetohydrodynamic wind, 
which is generated by the central pulsar, terminates 
at a shock, with the nebula identified as a postshock 
flow shining in synchrotron radiation. The central 
cavity of the nebula is occupied by the unseen wind. 
The shock is supposed to occur at the standing inner 
wisp. 

The KC model is very successful explaining the 
synchrotron luminosity, spectrum and frequency- 
dependent size of the nebula. An important conclu- 
sion of the KC model is that the energy of the wind is 
conveyed mostly by kinetic energy in the bulk motion 
of the plasma. Because the energy flux is in the form 
of an electromagnetic field at the base of the wind, 
this means that the efficiency of the wind acceleration 
is extremely high; KC found it to be 99.7%. 

The principal parameters of the pulsar wind are 
(1) its luminosity L^, (2) the Lorentz factor of 



the bulk flow and (3) the ratio a of the electromag- 
netic energy flux to the kinetic energy flux, which 
is referred to as the magnetization parameter. 
is essentially the spin-down luminosity « 5 x 10'^^ 
erg s~^. The remaining two parameters. and 
(T, together with the nebula pressure P/v, or equiv- 
alently the equipartition field Beg — V ^-'kPj^ , deter- 
mine the overall synchrotron spectrum. Conversely, 
the synchrotron spectrum tells us about the parame- 
ters. Given the synchrotron luminosity of 2 x 10"^^crg 
s~^, the nebula size of ~ 1 pc, and the peak and turn- 
off synchrotron spectrum energies of 2 eV and 10^ eV, 
respectively, one finds 7^, = 3.3 x 10^, cr = 3.8 x 10""^, 
and Bf.q = 0.38 mG. This result can be obtained 
even with an order-of-magnitude estimate (Shibata, 
Kawai and Tamura 1998). More rigorous fitting of 
the observed spectrum of the whole nebula gave sim- 
ilar values (e.g., KG, Atoyan and Aharonian 1996). 
The field strength has been confirmed by observa- 
tions of inverse Gompton emission in the TeV band 
(Wcckcs et al. 1989; Hillas ct al. 1998). Thus, the 
dominance of the kinetic energy of the wind seems 
very firm. 

From a theoretical point of view, however, the 
smallness of tr, or in other words dominance of the 
kinetic energy, is a mystery. No wind theory has 
been able to explain how such a high efficiency of 
acceleration is achieved. 

Ghandra observation clearly shows a disk-jet struc- 
ture and moving wisps with seeds of ~ 0.45c (Mori 
2002), where c is the speed of light. Because the KG 
model is spherically symmetric and steady, it may 
seem insufficient to understand the highly structured 
and dynamical nature seen by Ghandra. However, 
the basic idea that a kinetic dominant wind shocks 
and shines seems still firm and convincing. One may 
assume that the equatorial wind has different param- 
eters than the polar wind. Such a latitude depen- 
dence of the wind parameters may suffice to explain 
the apparent disk-jet structure although how such a 
latitude dependence is made is not known. 

In this paper, we suggest that high spatial reso- 
lution of Ghandra affords a chance to examine the 
assumptions which were made in the KG model but 
have yet to be checked. Among the assumptions, the 
ideal-MHD condition (no dissipation) and toroidal 



field approximation are of particular importance. If 
these assumptions are not adequate, the past conclu- 
sion of small cr may need to be reconsidered. 

We model the nebula in 3-dimensions based on the 
KG picture and reproduce an image, which can be 
compared with the Ghandra observation (Mori 2002) 
is made. We shall show that a considerable change 
to the KG picture is required to reproduce the Ghan- 
dra image. In this paper, we suggest that disordered 
magnetic field in the nebula is needed . Some pro- 
cess which converts magnetic energy into thermal and 
kinetic energy, such as magnetic reconnection, may 
take place in the nebula. In a subsequent paper, 
spatially-resolved X-ray spectra will be described and 
compared with Ghandra results. 

2 A 3D Model 

We postulate that the nebula flow obeys the KG 
steady solution and do not solve the dynamics. The 
properties of the KG flow are summarized as follows. 
If a is much less than unity as was suggested, the 
speed of the flow is ^ (l/3)c just after the shock 
and decreases rapidly with distance from the pulsar 
R as V (X R~^\ because the flow is subsonic, the 
pressure and density n are roughly uniform as for 
adiabatic expansion such that the mass conservation, 
nR^V wconst. (implying a decrease of the flow veloc- 
ity). Due to deceleration, the magnetic field accu- 
mulates and is amplified according to the frozcn-in 
condition, B (x r. Once the magnetic field grows as 
large as the equipartition field, the magnetic pressure 
becomes important in the flow dynamics. As a result, 
the flow speed saturates. This takes place where the 
nebula is brightest (at ~ shock radii). The 

smaller cr, the larger and brighter the nebula. The 
smallness of cr is then required to explain the lumi- 
nosity and the extent of the nebula. The indicated 
flow speed is small if cr is small: V/c ~ 3cr. It is no- 
table here that the above flow dynamics depends on 
the assumption of the ideal-MHD condition. 

The KG model is spherically symmetric and ob- 
viously inadequate to account for the observed mor- 
phology. We therefore restricted ourselves to an equa- 
torial region of the KG spherical model with half 
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Figure 1: The three-dimensional structure we as- 
sumed for reproduction of image. The spherical flow 
by Kennel and Coroniti (1984) is picked up for the 
disk and polar flows. 



width of Qeq ~ 10° and cut the intermediate latitude 
regions out, such that the disk may be reproduced in 
an image (see Fig. 1). Although we assume a polar 
flow by leaving the polar region with a semi-opening 
angle of Qpoi ^ 10° for reproduction of the jet image, 
this is just in an artist's spirit, and we do not provide 
any analysis of the polar jets in this paper. 

Our kinematic scheme to reproduce images of the 
nebula is made so that different types of the flow dy- 
namics can be applied in the future; i.e., for a given 
velocity field, we trace a fluid element and associ- 
ated particle distribution function according to the 
Lagrangian view of fluid. 

Given the radial velocity held V{t, R) as a function 
of time and radial distance from the pulsar, the posi- 
tion of the fluid element, the toroidal magnetic fleld, 
and the proper density are respectively obtained from 
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where t is the observer's time and F = (1— y^/c^)""'^/^ 



so uti' 



the Lorentz factor of the flow. We use the KG 
ion, which is given analytically, for the veloc- 
fleld. Although the KG solution does not include 
y|ichrotron losses, its effect on V is supposed to be 
all because synchrotron losses are about 10% of 
; kinetic energy of the flow. 

For the energy distribution of the particles, we 
agiin invoke the KG picture: a power law distribu- 
n is built up immediately after the shock, and the 
accelerated particles simply lose their energy 
the postshock flow by adiabatic and synchrotron 
looses. We trace the energy e (normalized by mc?) of 
;h particle in a fluid frame by 



h 3ck- 



(Ine) = (Inn) + - 

Dt'^ ' 3Dt'^ ' e 



loss 



where 



loss 



= -CTTCe C/mag, 



(4) 



(5) 



and the magnetic energy density, C/mag = B^/SttF^, is 
measured in the proper frame, and ctt is the Thomson 
cross section. The proper time t' is related to t by 
Dt/Dt' = F. 

The distribution function is deflned by 



dn = /(i, R\ e, 0) sinOdddipde, 



(6) 



where 6 is the pitch angle with respect to the local 
fleld and ip is the azimuth. Note that the distribution 
function is deflned in the flow proper frame. We as- 
sume that the postshock distribution follows a power 
law with index p (> 1) in between the minimum en- 
ergy emin and the maximum energy Cmax, and it is 
isotropic such that 



Me,) = ^n,e7\ 



(7) 



where the suffix 'i' indicates 'injection' at the post- 
shock region, and the normalization is given by 



(tmin/^max)^ 
In (^max/^min) 



p^i 
1 



(8) 



and n\ = JJJ f i{e) sin Bdddipde gives the postshock 
proper density, emax is assumed to be a maximum 
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attainable value, ei?2^s/"T.c^, where Rs is the shock 
distance from the pulsar, and B2 is the postshock 
field. Emin is determined so that the pressure calcu- 
lated from /i satisfies the shock jump condition. 

We solve Q numerically for a sample of particles 
in a given fluid element, and thereby we obtain e as 
a function of ei and t. Then we calculate distribution 
functions from 



/(e(i?)) = -MeO^. 



(9) 



For steady state models, integration for a single fluid 
element gives distribution functions in the whole neb- 
ula. 

It is obvious that the above kinematic scheme can 
be easily generalized for non-steady and non-radial 
flow, which can be obtained by numerical MHD sim- 
ulations. 



3 Reproduction of the Nebula 
Image 

3.1 Synchrotron Specific Emissivity 

Once the evolution equations Q-Q) are solved, and 
the distribution function is obtained by lO, it is 
straight forward to get volume emissivity, which is 
integrated to give a nebula image. 

Since the nebula flow is relativistic, a Lorentz 
transformation is applied between the flow proper 
frame and the observer's frame (or rather the pul- 
sar frame, in which the pulsar is at rest). Let us 
denote the 4- vector of a synchrotron photon by (lu / c, 
k) in the observer's frame and {uj' /c, k') in the proper 
frame. If the ideal-MHD condition E + V x B/c = 
holds in the nebula flow, the transformation of the 
electromagnetic field is simpler: 



E' = 

B\l = 
B', = 



0, (10) 
B\\ = {V ■ B)V/V^, (11) 
Bjr=iB-Bn)/T, (12) 



where the primes indicate the quantities in the flow 
frame and 11 and _L are based on the directed of the 



flow velocity V. There is no electric flcld in the 
plasma flow frame. 

The spectral power of a relativistic particle with 
pitch angle 9 (the angle of the particle motion to the 
local magnetic field in the proper frame) is 

Vsi{i^',e,9) ^ 2aTcUru.ge^smH Siuj';uj,). (13) 

For the monochromatic approximation, we use 
S{uj';lOc) — 5{lo' — cjc), and for the relativistic ap- 
proximation, 



^( , ^ 9^/3 /u;' 
S{uj ;a;c) = F — 



(14) 



r-OQ 

where F(x) = x Ks{^)d^, and 

J X 



3e|B'|e2 gin 6* 
2mc 



(15) 



is the critical frequency. The synchrotron power of 
a single particle is strongly beamed within a width 
of ~ e~^. Therefore, the emission into a frequency 
interval duj' and in a solid angle dfl' directed toward 
the observer is given by 

/>oo 

ju;'{9)duj'dn' / Vsi{uj',9,e)f{e,9)d€dn'duj', 
Jo 

(16) 

where 9 is given by cos 9 = n' ■ B' /\B'\, and n' is the 
unit vector directing to the observer in the proper 
frame. Below, n indicates the observer's direction in 
the observer's frame. 

For the link between the proper frame and the ob- 
server's frame, we include the Doppler effects, 

^' , f^' + P 
UJ = — -, ana 11= , 17 

r(i-/3^) ^ i + /3/i'' ^ ' 

where ^ — n-V and ji' = n' -V . The unit vector of 
the flow direction is denoted by V . The transforma- 
tion between the received power dPr and the emitted 
power dP' (Rybicki & Lightman 1979) is given by 



dP dP' 

-r3(i + /3/i')3^^- 



dP' 



dilduj 



dVl'duj' T'Hl^ I3^if dn'doj' 

(18) 
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where the Doppler effect (|17|l has been taken into 
account. FinaUy the Lorentz contraction is dN = 
rf{e,9)dfl'de, where N is the number density in the 
observer's frame. Thus the emissivity in the ob- 
server's frame becomes 



where 



jUn)^C / rsiiu;',e,e)f{e,6)de, (19) 



(20) 



3.2 Viewing Angle 

In order to specify the viewing angle of the observer, 
we relate the 'observer's coordinate' X = (X, Y, Z), 
where +X directed toward the observer and +Z di- 
rected toward north on the sky, to the 'nebula coor- 
dinate' X — {x, y, z), where the z-axis coincides with 
the symmetry axis of the nebula which is believed to 
be the rotation axis of the pulsar. We use 48° and 
28° as the position angle and the inclination angle of 
the symmetry axis, respectively. 

An image of the nebula is obtained by 



jUX,Y,Z,n)dX. 



(21) 



3.3 Summary of the procedure 



For a given observation frequency, the integration 
(|21|l is done numerically for each 'pixel' at {Y, Z). 
The integrand is calculated as follows: 

1. Given X = {X, Y, Z) and tj, 

2. transformation from the observer's coordinates 
to the nebula coordinates is done; for the po- 
sition, (X, y, Z) —f (x,y,z), and also for 
the componets of the observer's direction, n = 
(1,0,0) ^ n = {n^^,ny,n^). 

3. The flow velocity V{R) and the magnetic field 
B{R) at the point are obtained by using the KG 
solution. Note that the velocity and the mag- 
netic field are respectively radial and azimuthal 
in the 'nebula coordinate'. 



4. The observation frequency, the direction and the 
local magnetic field are transformed into those in 



5. The pitch angle of the particles directed toward 
the observer is obtained. 

6. Regarding the distribution function at (.t, y, z), 
the emissivity in the flow frame is calculated (we 
use the monochromatic approximation). 

7. Finally, the emissivity is converted to the volume 
emissivity juj{n) at {X, Y, Z) in the observer's 
frame by multiplying by the Doppler factor C. 

4 Results 

One may expect a ring-like structure in the repro- 
duced image, such as observed with Chandra (see the 
top panel of Fig. 2), since we have assumed that the 
flow is restricted within a disk. The expected radius 
of the ring will be ^ I/S^/ct (about 6 shock radii 
for a — 0.003), at which point the nebula brightens 
due to the amplified magnetic field. However, what 
we have is not a ring-like but is rather a 'lip-shaped' 
image shown in the bottom panel of Fig. 2. At the 
north-east and south-west corners of the expected el- 
lipse (ring), the pitch angles of the particles directed 
toward us are small, and therefore the surface bright- 
ness is reduced. This effect combined with the cen- 
tral cavity yields an image which is 'lip-shaped'. The 
smallness of pitch angle actually has two effects. One 
is that the single-particle emissivity is proportional to 
sin^ 9, which is small. The other is due to number of 
contributing particles. For a given observation fre- 
quency, the energy of the particles radiating at the 
frequency is higher for smaller pitch angles, so that 
the number of particles contributing to the frequency 
is smaller because of the negative slope of the distri- 
bution function. 

Another important point to consider is the inten- 
sity ratio between the front and back sides of the ring. 
We obtain a value of 1.3, but observed value is ~ 5 
(Pelhng et al. 1987, Willingale et al, 2001). The 
weak contrast is caused by deceleration of the neb- 
ula flow (i.e. the smallness of a). Mori et al. (2003) 
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Figure 3: Surface brightness along the semi-major 
axis of the torus in the north-east direction. The 
upper panel shows the calculation with different cr's 
while the lower panel is measured from the Chandra 
observation (Mori 2002). The distance is in units 
of the shock radius (which is supposed to be about 
13" — 14") for the model, and in arc seconds for the 
observation. The brightness is in units of 10^ erg s~-^ 
cm~^ str~^ eV~^ for the model and in counts s""'^ 
arcsec"^ for the observation. 



Figure 2: The Chandra image (top) and a reproduced 
image (bottom), where we assume a postshock flow 
with cr = 0.003 by Kennel and Coroniti (1984) but the 
flow is assumed to be restricted within an equatorial 
region. For the bottom image, the gray scale is in 
units of 0.016 erg s~^ cm~^ str eV. 



suggest that the ratio is about 3 with Chandra. This 
value is still incompatible with the KC picture. As 
long as the intensity contrast is attributed to Doppler 
boosting, such a weak contrast is unavoidable in the 
frame work of the KC model. 

How surface brightness changes with distance from 
the shock depends on a. Mori (2002) also measured 
the surface brightness along the semi-major axis of 
the torus, from which Doppler boosting should not 
affect the brightness. This result (bottom panel of 
Fig. 3) is compared with the present model (top panel 
of Fig. 3), for which we provide curves of varying 
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cr. The present model does not reproduce the first 
peak in the observation, which corresponds to the 
inner ring. However, it is notable that the brightness 
distribution of the inner ring is similar to that for the 
a = 0.1 model. The location of synchrotron burn off 
is reproduced by the a = 0.01 model. Finally, we 
point out that the surface brightness decreases much 
faster with distance in the observation than in the 
model. 

As indicated by the 'lip-shaped' image, the abso- 
lute value of the surface brightness is much less than 
observed along the semi-major axis. Because the re- 
produced image includes only the disk component, 
to which we restricted ourselves (rather than assum- 
ing the spherical KC model), the X-ray luminosity of 
the reproduced image is also smaller than the obser- 
vation. For the image in Fig. 2, we use parameters 
given by KC: = 5 x lO^^ erg a~^, = 3 x 10^^ 
cm, 7u, = 3 X 10^, and p = 3. In this case, we have 
I'Li, ^ 10^^ erg sec~^ at 1 keV. 

5 Discussion 

Applying the KC model, we reconstruct an X-ray im- 
age which is found to be inconsistent with the Chan- 
dra image. Owing to a pure toroidal field and uniform 
pitch angle distribution, the reproduced image is not 
ring-like but 'lip-shaped'. Furthermore, the surface 
brightness contrast between the front and back sides 
of the ring is much less than the observed. The 
weak contrast is simply due to the smallness of cr, 
by which the postshock flow slows down quickly after 
the shock. The assumptions of the toroidal field and 
the smallness of cr are thus found to be incompatible 
with the observation. 

If we assume isotropic emission in the proper frame 
such as is expected in a turbulent field, then the 
ring-like structure is reproduced as shown in Fig. 4. 
As would be expected, we find that such a turbu- 
lent component must be at least comparable to the 
mean toroidal field in order to reproduce the ring im- 
age. Although an another solution can be to adopt 
a contrived pitch angle distribution, we think this is 
unlikely. The image in Fig. 4 is produced in the fol- 
lowing way: (1) assume the magnetic field is random 



10 ^ 2>!lO ^ 3x10 ^ 4x10 ^ 

Figure 4: An image reproduced with assumptions of 
a turbulent field and a high speed fiow. See text in 
detail. 

so that the emissivity is isotropic in the proper frame; 
(2) set the fiow velocity to be 0.2c by hand, ignoring 
the fiow dynamics; (3) let the distribution function 
and the field strength follow the KC model. Thus, 
the random field and the fast flow are essential to 
reproduce the image. 

With this practice, we suggest that the nebula field 
is far from pure toroidal, but rather is disordered. 
Such a disordered field can be produced by magnetic 
reconnection or some instability of the toroidal fiux 
tubes. If there is dissipation of the magnetic field, 
the fiow dynamics is importantly changed, as is the 
fiow speed. Even if the value of cr, which is defined in 
the wind, is not small, dissipation in the nebula fiow 
causes deceleration and brightening. In this scncc, cr 
is effectively small so that the luminosity of the neb- 
ula will be explained as the KC model. But. such a 
simple heating may not always be good for explaining 
the surface brightness contrast because of the decel- 
eration. Recently, Komissarov & Lyubarsky (2003) 
provide an MHD simulation for the Crab Nebula, sug- 
gesting a complicated flow pattern and a high speed 
flow such that the brightness contrast can be repro- 
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duced. Three dimensional motions associated with 
magnetic energy conversion in the nebula will consid- 
erably change the picture of the nebula. Numerical 
simulations for the nebular flow is of particluar im- 
portance in the future study. As noted, the present 
image-production scheme will be easily extended to 
combine with such mumerical results. 

A model explaining the Chandra observation may 
be constructed if we assume a larger a and a sub- 
sequent magnetic energy conversion into heat and 
plasma kinetic energy, such as magnetic reconnec- 
tion, in the nebula flow. Suppose a is rather large, 
then the postshock flow must be faster. The inner 
ring is formed at the shock. The brightness distribu- 
tion will be similar to that of cr = 0.1 in Fig. 3. As the 
flow proceeds outward, the magnetic energy conver- 
sion takes place (accelerating and heating the flow). 
This causes the second brightening. Subsequently, 
the synchrotron burn-off provides the outer boundary 
of the torus. The Doppler effect will cause a higher 
brightness contrast. We note that the smallness of a 
is not obvious if non-ideal-MHD is introduced in the 
nebula flow. 

The above hypothesis explains the disk formation. 
Obliqueness of the pulsar causes a series of current 
sheets with an interval of the light cylinder radius 
lO^cm) in the equatorial region. If the current 
sheets dissipate in the nebula, the synchrotron emis- 
sion brightening is restricted in the equatorial region, 
where reconnection takes place. The appearance of 
pulsar nebulae should depend on obliqueness of in- 
dividual pulsars. High obliqueness results in a thick 
disk and a high efflciency of synchrotron luminosity, 
while near-alignment causes a faint nebulae. 

The possibility of a dissipative process in the neb- 
ula may be examined more rigorously with spatially- 
resolved spectra, for which we will compare the model 
and the observation in a subsequent paper. 
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